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ABSTRACT 

We propose a unified model to explain Quasi-Periodic-Oscillation (QPO), particularly of high fre- 
quency, observed from black hole and neutron star systems globally. We consider accreting systems to 
be damped harmonic oscillators exhibiting epicyclic oscillations with higher order nonlinear resonance 
to explain QPO. The resonance is expected to be driven by the disturbance from the compact object 
at its spin frequency. The model explains various properties parallelly for both types of the compact 
object. It describes QPOs successfully for ten different compact sources. Based on it, we predict the spin 
frequency of the neutron star Sco X-l and specific angular momentum of black holes GRO J1655-40, 
XTE J1550-564, H1743-322, GRS 1915+105. 

Subject headings: black hole physics — stars: neutron — gravitation — relativity 



1. INTRODUCTION 

The origin of Quasi-Periodic-Oscillation (QPO) and its 
properties are ill-understood phenomena in modern astro- 
physics. Observed QPO frequencies from black hole, neu- 
tron star, as well as white dwarf systems range from milli 
Hertz (mHz) to kilo Hertz (kHz) with varying properties. 
Among them most interesting and puzzling QPOs are the 
kHz ones. To date, more than 20 accreting neutron stars 
in Low Mass X-Ray Binaries (LMXB) have been found 
exhibiting such kHz frequency QPOs (van der Klis 2000; 
2006). Several black holes have also been found exhibiting 
the high frequency (HF) QPOs of a few hundreds of Hz. 
The most puzzling cases are those when a pair of QPO 
forms, particularly in neutron star systems, and QPOs in 
the pair appear to be separated in frequency either by the 
order of the spin frequency of the neutron star, apparently 
for slow rotators, or by half of the spin frequency, for fast 
rotators. An example of the later cases is 4U 1636-53 with 
the spin frequency v s ~ 580 Hz, while of the former ones is 
4U 1702-429 with v s ~ 330 Hz. Moreover, in several occa- 
sions the frequency separation decreases with the increase 
of one of the QPO frequencies. In black hole binaries also 
the HF QPOs have been found in pairs, but seem to appear 
at a 3 : 2 ratio (Remillard et al. 2002) that favors the idea 
of a resonance mechanism behind their origin. The exam- 
ples are GRO J1655-40, XTE J1550-564, GRS 1915+105 
which exhibit the HF QPOs. The difference of the twin 
HF/kHz (van der Klis 2005) QPOs between black holes 
and neutron stars might be due to presence of a rotating 
magnetosphere in the later cases imprinting directly the 
spin frequency on the oscillations of the disk. However, 
sometimes QPOs from neutron stars, e.g. Sco X-l, also 
seem to appear at a 3 : 2 ratio (Abramowicz et al. 2003). 

Magnetically driven warping and precession instability 
are expected to occur in the inner region of an accretion 
disk (Lai 1999) by interaction between the surface current 
density in the disk and the stellar magnetic field. This 
produces motions at low frequencies, much lower than 
the orbital frequency accounting for the low frequency 



QPOs. This is similar to the Lense-Thirring precession 
shown by Stella & Vietri (1998). The mHz QPOs in some 
systems can be explained with magnetic warping process 
for strongly magnetized neutron stars (Shirakawa & Lai 
2002). Titarchuk (2003), on the other hand, described 
QPO by the Rayleigh- Taylor instability associated with 
Rossby waves and rotational splitting. The effect of non- 
linear coupling between g-mode oscillations in a thin rela- 
tivistic disk and warp was examined by Kato (2004) for a 
static compact object. 

A strong correlation between low and high frequency 
QPOs has been found by recent observations (Psaltis, Bcl- 
loni & van der Klis 1999; Mauche 2002; Belloni, Psaltis & 
van der Klis 2002) of accretion disks around compact ob- 
jects. Holding the relation over mHz to kHz range strongly 
supports the idea that QPOs are universal physical pro- 
cesses, independent of the nature of the compact object. 
Observing universal nature between the HF/kHz QPOs in 
black holes and neutron stars, McClintock & Remillard 
(2004) suggested that they are due to the orbital motion 
of the accreting matter near the marginally stable circu- 
lar orbit. Titarchuk & Wood (2002) explained the cor- 
relation in terms of a centrifugal barrier model (see also 
Titarchuk, Lapidus & Muslimov 1998). Indeed, earlier by 
numerical simulation, Molteni, Sponholz & Chakrabarti 
(1996) showed QPO in the black hole accretion disks from 
the shock oscillation induced by the strong centrifugal bar- 
rier. Later, Chakrabarti & Manickam (2000) verified this 
analyzing observational data. 

In early, it has been shown that QPOs may arise from 
nonlinear resonance phenomena in an accretion disk placed 
in strong gravity (e.g. Kluzniak et al. 2004; Petri 2005; 
Blacs et al. 2007). This is particularly suggestive because 
of their variation in frequency, while we know that the 
varying frequency in a characteristic range is a well studied 
feature of nonlinear oscillators. Indeed accretion dynam- 
ics and the disk structure are determined by the nonlin- 
ear hydrodynamic/magnetohydrodynamic equations with 
the effective gravitational potential which is not harmonic. 
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Therefore, a nonlinear response is expected. The above au- 
thors argued for the resonance between epicyclic motions 
of accreting matter around a rotating compact object re- 
sulting QPOs, particularly the HF/kHz ones. However, 
the separation between vertical and radial epicyclic fre- 
quencies as a function of radial coordinate increases mostly 
with the increase of cither of the frequencies. Although 
the separation of QPO frequencies in a pair observed from 
XTE J1807-294 appears constant (Linares et al. 2005) and 
from Aql X-l seems increasing (Barret, Boutelier & Miller 
2008) with the increase of either of the frequencies, in most 
of the cases it is opposite (modeled by Stella & Vietri 
1999; Karas 1999) which contradicts the model based on 
the resonance between epicyclic motions. At radii smaller 
than the radius showing peak of radial epicyclic frequency, 
while the separation decreases with the increase of radial 
epicyclic frequency, the frequency separation is too large, 
much larger than the spin frequency of the star. Indeed, 
based on Cir X-l data, Belloni, Mendez & Homan (2007) 
argued against the model. Therefore, one has to invoke 
new physics to describe multiple QPO features. 

In the present paper, we propose a global resonance 
model based on higher order nonlinearity to describe black 
hole and neutron star QPOs together. The model success- 
fully describes important features including variation of 
the QPO frequency separation in a pair as a function of 
QPO frequency itself. On the mission of reproducing ob- 
served QPOs, the model predicts the spin parameter (spe- 
cific angular momentum or the Kerr parameter) of black 
holes as well as the spin frequency of a neutron star. The 
paper is organized as follows. In the next section, we es- 
tablish the model describing the underlying physics briefly. 
In §3 and §4 we reproduce various properties of observed 
kHz/HF QPOs from six neutron stars and four black holes 
respectively, minimizing free parameters. Finally, we end 
with a summary in §5. 
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(apart from all the primary and secondary resonance con- 
ditions cv* = dv m , where v m corresponds to the mth 
mode, a, bi, c, d are integers) such that 



a + ^l&il =j, 



(2) 



i=i 



where j = k + 1 when k is the order of nonlinearity. Cer- 
tainly the order of nonlinearity determines the type of res- 
onance conditions. For details, see Nayfeh & Mook (1979). 

Now we consider an accretion disk to be a system of 
higher degree nonlinearity and described as a damped os- 
cillator with higher degree anharmonicity. Then we con- 
sider the possible resonance therein. The resonance is 
expected to be driven by the combination of the strong 
disturbance by the compact object with spin v s and other 
already existent (weaker) disturbances in the disk particles 
(and hence disk) at the frequencies of the radial (v r ) and 
vertical (y z ) epicyclic oscillations given by (e.g. Shapiro & 
Teukolsky 1983) 



v r = y^A^VF-a) 2 , v z = ^^Jr 2 - las/i + 3a 2 (3) 



where A 



2r + a 2 , r is radial coordinate expressed in 



unit of GM/c 2 , a is specific angular momentum (spin pa- 
rameter) of the compact object expressed in unit of GM/c, 
G the Newton's gravitation constant, c is speed of light, M 
is mass of the compact object, and v Q is orbital frequency 
of the disk particles given by 



2itv = Q, 



r 3 / 2 + a GM ' 



(4) 



2. MODEL 

Wagoner (1999) and Kato (2004) investigated small os- 
cillations of geometrically thin accretion disks around a 
rotating compact object in the linear regime in the con- 
text of QPO. Later, Rezzolla et al. (2003) studied the 
radial oscillations of twodimensional accretion tori in the 
Schwarzschild spacetime. Kluzniak et al. (2004) investi- 
gated the response to a transient external perturbation of 
an ideal gas torus around a nonrotating black hole or neu- 
tron star. They found that the torus performs harmonic 
global oscillations in both the radial and the vertical di- 
rections with a variable amplitude in the later direction. 
These effects indicate the existence of a nonlinear coupling 
into the system. 

We know that a system of N degrees of freedom has 

N linear natural frequencies denoting v\,vi,v-z, 

(Nayfeh & Mook 1979). Depending on the order of non- 
linearity to the system, these frequencies are commensu- 
rable with certain relationships which may cause the cor- 
responding modes to be strongly coupled and yield an in- 
ternal resonance. If the system is excited by a mode of 
frequency v*, then frequencies might have the commensu- 



We describe the system schematically in Fig. 1. It can 
be seen to have composed of two oscillators, one in the 
radial direction (radial epicyclic motion) and other in the 
vertical direction (vertical epicyclic motion), with differ- 
ent spring constants. Note that spring constants here vary 
with radial coordinate. The basic idea is that the mode 
corresponding to v s arised due to the disturbance created 
by rotation of the compact object will couple to the ones 
corresponding to v T and v z , already exist in the disk, excit- 
ing new modes with frequencies v r ^ z ±pv s (Nayfeh & Mook 
1979), where p is a number, e.g. n/2 when n is an integer 
(1, 2 for the present model). Now at a certain radius in the 
nonlinear regime where v s /2 (or v s in the linear regime) 
is close to v z — v r (Nayfeh & Mook 1979; Kluzniak et al. 
2004), if v s /2 (or v s ) is also coincidentally close to the fre- 
quency difference of any two newly excited modes, then a 
resonance may occur which locks the frequency difference 
of two excited modes at v s j2 (or v s ). 

For a neutron star having a magnetosphere coupled with 
the disk matter, the mode with v s can easily disturb disk 
particles through the interaction between surface current 
density in the disk and the stellar magnetic field and get 
coupled with the modes existing in the disk. Ghosh & 
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Lamb (1979a, b) showed, close to the stellar surface, how 
does the magnetic field channel matter to the polar caps. 
Later, Anzer & Borner (1980, 1983) showed the influence 
of the Kelvin-Hclniholtz instability. However, LMXB have 
a very low surface field strength and a low critical fastness 
parameter compared to highly magnetized neutron stars 
(Li & Wickramasinghe 1997; Li & Wang 1999). Therefore, 
in this case, the strength of the disturbance is expected to 
be low and the accretion disk usually extends close to the 
stellar surface (Popham & Sunyaev 2001). 

For a black hole, on the other hand, energy and an- 
gular momentum can be transferred with the magnetic 
field, associated with its disk, threading the horizon and 
connecting to the surrounding environment which may be 
looked as a variant of Blandford-Znajek process (Bland- 
ford & Znajek 1977). In fact, this kind of magnetic cou- 
pling process was already supported by XMM-Newton ob- 
servations (e.g. Wilms 2001; Li 2002). This confirms the 
possible excitation of new modes in the disk around a spin- 
ning a black hole. 

Note, as remarked by Karas (1999) in a similar context, 
that the disk particle size should be finite to visibly modu- 
late the observed radiation due to the resonance for a long 
period. However, shearing effect destroys large particles. 
Abramowicz et al. (1992) suggested that formation of vor- 
tices in accretion disks could help in surviving finite size 
particles which was verified by Adams & Watkins (1995). 
While the origin of vortex is not fully understood yet, finite 
amplitude threcdimensional secondary perturbations may 
form elliptical vortices in disks surviving for a long period 
(Mukhopadhyay, Afshordi, Narayan 2005; Mukhopadhyay 
2006). In fact, Mukhopadhyay (2008), based on it, argued 
for the origin of turbulent viscosity in accretion disks. 

Therefore, we rewrite (1) for an accretion disk as 

(n - m+ l)v s = b\v r + b 2 v z (5) 

which leads to 

, n v s m 

V r + ^V a = V s - — + jVs (6) 

with — bi = 62 = 2 and m, n are integers. Now we propose 
the higher and lower QPO frequencies of a pair respec- 
tively to be 

i n v s 
Vh = v r + -v s , vi = v z -—. (7) 

Hence, from (2) and (6) we understand the possible or- 
der of nonlinearity in accretion disks exhibiting QPOs 
k = n — m + 4. In the next section, we will reproduce 
several observational results, were not addressed by other 
QPO models, what justify the proposition. 

In the disk around a neutron star, at an appropriate 
radius where v z — v r ~ u s /2 (Kluzniak et al. 2004) and 
the resonance is supposed to take place, for n = m = 1, 
Av = Vh — ui ~ v s /2, and for n = m = 2, Av ~ v s . 
n = 1 corresponds to a nonlinear coupling between the ra- 
dial epicyclic mode and the disturbance due to spin of the 
neutron star, which results in Av locking in the nonlinear 
regime with m = 1. On the other hand, n = 2 corre- 
sponds to a linear coupling, which results in Av locking 

2 Note that the strength of the disturbance is expected to be higher 
hence, plausibly a nonlinear resonance is in the former case and a linea 
relation between QPO and spin frequencies of neutron stars support th 



in the linear regime with m = 2. Hence, for the reso- 
nance n must be equal to m. The resonance with, e.g., 
n = 2, m = 1 for a neutron star corresponds to a linear 
coupling of the epicyclic and the disturbance modes. This 
results in the locking of Av (~ v s /2) to be in the nonlinear 
regime which is forbidden or very weak to observe. Sim- 
ilarly, n = 3, m = 3 may correspond to higher harmonics 
of the coupling in the nonlinear regime which are again 
expected to be weak to make any observable effects. It is 
expected that for a fast rotating neutron star the distur- 
bance occurs in the nonlinear regime 2 . However, both v r 
and v z and then and v\ vary as functions of radial coor- 
dinate. Therefore, Av = Vf t — vi varies as well. We argue 
that in a small range of radius, as given in Table 1, the 
resonance condition remains satisfied approximately, when 
Av decreases slowly with v\ , hence explains observed data. 
Indeed observations show Av (may vary over 100Hz for a 
particular neutron star) to be of the order of v or v/2 only, 
but not to be v or v/2 exactly. In Fig. 2 we describe these 
facts for stars with Av ~ v s /2 and v s both. The vertical 
dot-dashed line indicates the location of the resonance. 

For a black hole, however, in absence of a magneto- 
sphere, the disturbance and then corresponding coupling 
between modes may not be nonlinear and occurs with the 
condition v z — v r £ v s , which results in the resonance lock- 
ing at the linear regime with n = m = 2 which pro- 
duces Av v s (sometimes <~ 2^ s /3) as shown in Figs. 
3a, b. However, if we enforce the resonance to occur at 
marginally stable circular orbit, at the very inner edge of 
an accretion disk where the underlying physics is expected 
to be nonlinear than that in the relatively outer edge, then 
a nonlinear resonance occurs with Av ^ v s /2 (sometime 
~ v s /5) and v z — v r £ v s for n = m = 1, as shown in Figs. 
3c, d. In the following sections this is discussed in detail. 

The difference between resonance conditions to form 
QPOs around black holes and neutron stars might be due 
to difference in their environment. For a neutron star, 
there is a rotating magnctosphcrc imprinting directly the 
spin frequency on the oscillations of the disk. However, en- 
ergy and angular momentum of a black hole with magnetic 
field connecting to the surrounding disk can be transferred 
to the disk through the variant Blandford-Znajek process. 
This leads to the strong resonance, expected to be driven 
by the disturbance at the spin frequency of the black hole, 
even though the exact physical mechanism remains to be 
determined. Another important point to note is that for 
a black hole there is no strict definition of v s correspond- 
ing to its boundary layer related to the spin parameter a, 
while for a neutron star it is obvious due to presence of its 
hard surface. Here we assume that the effective boundary 
of a black hole is at marginally stable circular orbit r ms 
and v s corresponds to the frequency of a test particle at 
r ms [see equations (9), (10)]. 

From observational data we can determine the spin fre- 
quency of a neutron star. In computing QPO frequencies 
from our model, we need to determine the spin parameter 
a. If we consider a neutron star to be spherical in shape 
with the equatorial radius R, spin frequency v s , mass M, 
radius of gyration Rq 1 then moment of inertia and the 

for a fast rotating neutron star compared to a slow rotating one and, 
• resonance in the later one. Indeed observations showing the apparent 
; argument. 
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spin parameter are computed 

I = MR%, a=!^r, (8) 



where ft s — 2-kv s . We know that for a solid sphere R 2 -, = 
2 i? 2 /5 and a for hollow sphere R G — 2 i? 2 /3. However, the 
shape of a very fast rotating neutron star is expected to be 
deviated from spherical to ellipsoidal. Moreover, neutron 
stars are not expected to be perfect solid bodies. Hence, in 
our calculation we choose 0.35 < (Rg/R) 2 < 0.5 in most 
cases (see e.g. Bejger & Haensel 2002; Cook, Shapiro & 
Teukolsky 1994). 

On the other hand, for a black hole of mass M, a is the 
most natural quantity what we supply as an input. Corre- 
sponding angular frequency of a test particle at any radius 
r in the spacetime around it is then given by (e.g. Shapiro 
& Teukolsky 1983) 

O _ 94>t _ 2 « /qx 

300 H +ra 2 + 2a 2 

We now define 

c 3 

O s = 2wu s = Q. B H{r = r ms ) (10) 

where r ms is marginally stable circular orbit and light in- 
side r ms is practically not expected to reach us. Therefore, 
supplying a we can determine v s . 

3. PROPERTIES OF NEUTRON STAR QPOS 

For at least seven of neutron stars exhibiting twin kHz 
QPOs we know the spin frequency (e.g. Yin et al. 2007; 
Boutloukos & Lamb 2008) In most occasions, QPO fre- 
quencies vary with time and show that the frequency dif- 
ference in the pair decreases with the increase of lower 
QPO frequency. Moreover, the frequency difference is of 
the order of half of the spin frequency of the star for fast 
rotators and of the spin frequency itself for slow rotators. 
We will show that our model successfully describes all the 
properties matching with observed data. In obtaining re- 
sults, we choose primarily M = 1.4M Q , (R G /R) 2 = 0.4 
along with a guess radius R; all are free parameters. 
However, this choice does not suffice observation mostly. 
Hence, subsequently, we discuss results with other values 
of parameters. 

3.1. Fast rotators 

We analyze QPOs from three fast rotating neutron stars 
whose spin frequencies are known: KS 1731-260 (Smith, 
Morgan & Bradt 1997), 4U 1636-53 (Jonker, Mendez, & 
van dcr Klis 2002) and 4U 1608-52 (Mendez et al. 1998). 
For various sets of input values of M and R G /R 2 we pre- 
dict radius of the neutron star, given in Table 1, in re- 
producing observed data by our theory, depicted in Figs. 
4a-c. KS 1731-260, till date, has shown to exhibit only 
one pair of QPO frequency established by our model eas- 
ily as depicted in Fig. 4a. However, the resulting R is 
large/unrealistic if mass of the star is considered to be 
M = 1.4M Q . If we choose the star not to be a solid sphere 
or to be an ellipsoid with (Rg/R) 2 = 0.5 or 0.35, then 
the resulting radius decreases to R = 16.4km which is in 



accordance with realistic equations of state (EOS) [Fried- 
man, Ipser & Parker 1986, hereafter eosl; Cook, Shapiro 
& Teukolsky 1994, hereafter eos2]. Figures 2a,b show the 
resonance to occur at rQPo =8 (indicated by the vertical 
dot-dashed line) for n = m = 1 with u z — v r ~ IS.V ~ v s j2. 

For 4U 1636-53 and 4U 1608-52, on the other hand, twin 
peaks are observed in several occasions. For 4U 1636-53, 
our theory has an excellent agreement with observation 
with x 2 ~ 28, as shown in Fig. 4b. If we exclude only data 
point falling outside the theoretical prediction, which does 
not appear to follow the same trend as other points follow 
and hence the significance of its detection appears to be 
less, then x 2 ~ 6. If either (Rg/R) 2 or M of the star is 
considered to be deviated from 0.4 or 1.4M Q respectively, 
then our theory reproduces observed data in accordance 
with realistic EOS (eosl, eos2). While at smaller z^-s, re- 
sults for 4U 1608-52 deviate from observed data, as shown 
in Fig. 4c, the predicted radii for all the parameter sets, 
except that with R — 10km, are well within the conven- 
tional range (eosl, eos2), given in Table 1. While the 
estimated x 2 for 4U 1608-52 is large (> 100) for the entire 
set of observed data, it is only ~ 8 when we fit a part of 
it, discarding remaining data points from the computation 
as of 4U 1636-53, by the dot-dashed line shown in Fig. 4c. 
However, the important point to note is that the number 
of data points in either of the cases is poor which may 
question the reliability of \ 2 values. Similar resonance as 
of KS 1731-260 also occurs for these sources as well. 

3.2. Slow rotators 

Results for slow rotating neutron stars whose spin 
frequencies are known: 4U 1702-429 (Markwardt, 
Strohmayer & Swank 1999), 4U 1728-34 (van Straaten et 
al. 2002; Mendez & van der Klis 1999), are depicted in 
Figs. 4d,e. For the sets of input values of M and (Rg/R) 2 
we predict various possible radii of the neutron star, given 
in Table 1. If we consider M = 1.4M , then the predicted 
R is too large with the choice of a spherical star. If the star 
is considered to be an ellipsoid and/or not to be a solid 
sphere and/or to have mass M < 1.4M Q , then R decreases 
to conventional values (eosl, eos2). On the other hand, all 
the parameter sets for 4U 1702-429, given in Table 1, cor- 
respond to high specific angular momentum of the neutron 
star. Therefore, we predict that the star 4U 1702-429 can 
not be a solid sphere and thus corresponding sets of pa- 
rameters in rows one and three appear to be ruled out. 
Figure 4d shows a perfect agreement of our theory with 
observation with x 2 < 1. The spin parameter for cases of 
4U 1728-34, given in Table 1, ranges 0.5 < a < 0.75. The 
parameter set with a <~ 0.75 given in the corresponding 
row one appears to be ruled out when the star is not ex- 
pected to be a solid sphere and the mass-radius combina- 
tion does not follow any realistic EOS. While the estimated 
X 2 for 4U 1728-34 is large (> 100) for the entire set of ob- 
served data, it is only ~ 7 when we fit a part of it, as of 
4U 1608-52, by the dot-dashed line shown in Fig. 4e. How- 
ever, as mentioned earlier, the poor number of data points 
may question the reliability of x 2 values. In Figs. 2c, d 
we describe the resonance and corresponding formation of 
QPOs for 4U 1702-429. Clearly, the resonance occurs in 
r QPO = 10 — 11, indicated by the vertical dot-dashed line, 
for n = m = 2 when v z — v r <~ Av/2 <~ v s /2. Similar 
resonance occurs for other slow rotators. 
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3.3. Estimating the spin frequency of Sco X-l 

The spin frequency of Sco X-l is not known yet. This 
source exhibits a pair of kHz QPO frequency with the sep- 
aration varies in a range ~ 225 — 310 Hz (Mendez & van 
der Klis 2000). We show in Fig. 4f the observed vari- 
ation of the frequency separation as a function of lower 
QPO frequency and compare it with that obtained from 
our model. Considering all possible sets of input parame- 
ters (along with v s which is now unknown) to obtain the 
best fit of observations, we present some of them in Ta- 
ble 1. We find that mass of Sco X-l must be less than 
IAMq and the sets of inputs with smaller v s and M fit 
the observed data better and argue that Sco X-l is a slow 
rotator with v s ~ 280 — 300. For the parameter set with 
v s = 292, x 2 ~ 66 when the number of observed data 
points is 39. However, if we recompute it discarding data 
points falling far outside the theoretical prediction, as of 
4U 1636-53, 4U 1608-52, 4U 1728-34, given by the dot- 
long-dashed curve, then \ 2 <~ 38 with observed data points 
35, assuring the excellence of the fitting. 

3.4. Quality factor 

Quality factor (in short Q-factor) of observed QPO fre- 
quencies has been found to be as large as 200. Any viable 
model for QPO should be able to reproduce this high Q- 
factor. Figure 5 describes Q-factor of the 775Hz QPO ob- 
served from 4U 1636-53, based on our model. The power 
depends mainly on mechanical resistance of the system 
along with cpicyclic frequencies 3 and frequency of the neu- 
tron star. Following the principles of a forced vibrating 
system of coupled oscillators we obtain the power approx- 
imately shown in Fig. 5a. However, Fig. 5b shows that the 
power falls off rapidly away from the resonance radius, be- 
cause the component of radial epicyclic frequency merges 
with the vertical epicyclic frequency with the increase of 
r. Therefore, the resonance condition remains satisfied 
roughly in a (narrow) range of radii, as given in Table 1. 
The change in amplitude of the power in that range of 
radii results in the observed variation of QPO frequency. 
Figure 5b shows that the QPO frequency of 4U 1636-53 
varies in the range tqpo = 7.6 — 8.5 (see Table 1) with 
amplitude of the power decreases upto one-third of peak 
amplitude. Similar features can also be found for other 
compact sources. 

The value of mechanical resistance r m of a system con- 
trols Q-factor. From the perturbed Navier-Stokes equa- 
tion, it is easy to understand that r m for an accretion 
disk, where QPOs occur, is proportional to the radial 
velocity gradient dv r /dr of the background flow. While 
for a Keplerian accretion disk (Shakura & Sunyaev 1973) 
dv r /dr is very small, for a sub-Keplerian advective disk 
(e.g. Mukhopadhyay 2003) it could be high which is dif- 
ficult to perturb rendering high mechanical resistance. In 
Fig. 5c we show the variation of Q-factor with r m . It 
clearly argues that a very large Q-factor is possible for 
QPOs from a Keplerian disk. 

4. PROPERTIES OF BLACK HOLE QPOS 

The pairs of observed HF QPO from black holes seem 
to appear at a 3 : 2 ratio. Mass (or possible mass range) 

3 Note that the spring constants of the underlying oscillators, which 
frequencies. 



of several such black holes, e.g. GRO J1655-40 (Orosz & 
Bailyn 1997; Shahbaz et al. 1999), XTE J1550-564 (Orosz 
et al. 2002), GRS 1915+105 (Greiner, Cuby & McCaugh- 
rcan 2001), H1743-322 (Miller et al. 2006) is already esti- 
mated from observed data. However, the spin of them is 
still not well established. We estimate the spin parameter 
of black holes in describing their observed QPOs by our 
model given in Table 2. 

We supply observed M (or range of M, whichever 
is available) and arbitrary values of a as inputs of our 
model to reproduce observed QPOs. Without knowledge 
of exact mass for GRO J1655-40, XTE J1550-564 and 
GRS 1915+105, the estimated range of mass from ob- 
served data is supplied for our computation and we find 
the corresponding suitable range of input parameter a 
reproducing observed QPOs. However, mass of H1743- 
322 appears to be less uncertain with an estimated value 
<~ 11M Q . So we use that value and find a suitable spin 
parameter. For clarity, we also choose a (hypothetical) 
range of mass around ~ 11M Q for this black hole and 
check in what extent a could vary. Our theory is able to 
produce QPOs at a 3 : 2 ratio with their observed values 
for n = m, = 2 at a radius outside r ms . Figures 3a, b de- 
scribe the resonance and corresponding formation of QPOs 
for the black hole GRO J 1655-40 with M = 6M and 
v s = 218.48Hz. The vertical dot-dashed line indicates the 
location of the resonance. It clearly says that v z — v r ,$ v s 
a t r QPO = 4.93 and the resonance occurs with Av ,$ v s 
(say, ~ 2v s /3, as shown in Fig. 3b). Similar resonance 
happens for other black holes. The corresponding spin pa- 
rameter ranges from ~ 0.6 to <~ 0.8, depending on mass of 
the black hole. 

However, if we enforce QPOs to happen at r ms strictly 
by changing mass (within observed range of estimate) and 
the spin parameter of a black hole, then they produce for 
n = m = 1 at a higher a. Figures 3c, d describe such 
a formation of QPOs for the source GRO J1655-40 with 
M = 7.05M Q and v s = 802.48Hz. It again implies that the 
resonance corresponds to v z — v r ^ v s , but with Av £ v s /2 
(say, <~ v s /5, as shown in Fig. 3d). However, in this case 
the required mass of black holes to reproduce observed 
QPOs may be out of the range estimated from observed 
data and the corresponding a appears to be similar/same 
for all four black holes. Interesting fact to note is that the 
cither case favors the fourth order nonlinearity in accre- 
tion disks around black holes and none of the parameter 
sets predict an extremally rotating black hole. 

In early, employing fully rclativistic accretion disk 
model, the spin of black holes was estimated based on 
a spectral analysis of the X-ray continuum (Shafee et 
al. 2006; McClintock 2006). While their estimated 
spin for GRO J1655-40 is in accordance with ours, for 
GRS 1915+105 it is higher (a > 0.98). However, Mid- 
dleton et al. (2006), based on a simple multicolor disk 
blackbody including full radiative transfer as well as rel- 
ativistic effects, estimated spin of GRS 1915+105 to be 
~ 0.7 which tallies with ours. Perhaps, different estimate 
is due to difference in method and the disk model used 
for computation. Indeed, the best model to describe a 
geometrically thick accretion disk is still a matter of con- 
troversy. In fact, Wang, Ye & Huang (2007) showed how 

vary with radial coordinate, arc related to the components of epicyclic 
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does the spin estimate vary with different method. 
5. SUMMARY 

We have described QPOs, particularly the kHz/HF 
ones, observed from several neutron stars and black holes 
(altogether ten compact sources), by a single model and 
predicted their spin. The model has addressed the vari- 
ation of QPO frequency separation in a pair as a func- 
tion of the QPO frequency itself of neutron stars, ob- 
served in several occasions. We argue that QPO is a result 
of higher (fourth) order nonlinear resonance in accretion 
disks. The model predicts the most possible set of pa- 
rameters: R ~ 15 — 16km, M ~ 1 — 1.2M Q for fast rota- 
tors; R ~ 16 - 17km, M ~ 0.8 - 1M for slow rotators, 
when (Rq/R) 2 ~ 0.4, that can explain all the observed 
QPO pairs. Discussing QPOs for five different neutron 
stars (KS 1731-260, 4U 1636-53, 4U 1608-52, 4U 1702-429, 
4U 1728-34) whose spin frequencies are known, we have 
predicted the spin frequency of Sco X-l which is unknown 
until this work. We argue that this is a slow rotator. 

We have addressed QPOs from four black 
holes (GRO J1655-40, XTE J1550-564, H1743-322, 
GRS 1915+105) whose mass (or range of mass) has already 
been predicted from observed data and QPO pair seems 
to appear at a 3 : 2 ratio. Based on the present model, 
we have predicted their spin parameters (o-s) which are 
not well established yet. According to the present model, 
none of them is an extremally rotating black hole. As 



our model explains QPOs observed from several compact 
objects including their specific properties, of that in black 
holes and neutron stars both, it favors the idea of QPOs 
to originate from a unique mechanism, independent of the 
nature of compact objects. A successful model should 
be able to reproduce properties globally as the one does 
presented here. 

Now the future job should be to perform numerical simu- 
lations to establish the resonance described here to confirm 
the model. Importantly, one should compare the energy 
released from the oscillation with that from rest of the disk 
(and the neutron star surface, if QPO is from a neutron 
star) to address visibility of oscillations. In this connec- 
tion, the variation of the disturbance in the disk due to 
rotation of the compact object should be employed prop- 
erly into the model to pinpoint the largest possible radius, 
beyond that QPO does not occur, accurately. 
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Table 1 

Physical parameters of neutron stars 



fast rotator 




M 


(Rg/RY 


n. m 


R 


range of tqpo 


KS 1731-260 


525.08 


1.4 


0.4 


1 


22.01 


one pair at 6.6 


KS 1731-260 


525.08 


1.1 


0.5 


1 


16.4 


one pair at 8 


KS 1731-260 


525.08 


1.1 


0.4 


1 


18.32 


one pair at 8 


KS 1731-260 


525.08 


0.9 


0.35 


1 


16.4 


one pair at 9.3 


4U 1636-53 


581.75 


1.4 


0.4 


1 


18.4 


6.6 - 7.4 


4U 1636-53 


581.75 


1.4 


0.49 


1 


16.8 


6.7 - 7.5 


4U 1636-53 


581.75 


1.2 


0.4 


1 


16.0 


7.5 - 8.4 


4U 1636-53 


581.75 


1.18 


0.35 


1 


16.8 


7.6 - 8.5 


4U 1608-52 


619 


1.4 


0.4 


1 


15.0 


7.0 - 8.9 


4U 1608-52 


619 


1.4 


0.5 


1 


13.5 


7.0 - 8.9 


4U 1608-52 


619 


1.1 


0.4 


1 


10.0 


8.5 - 10.7 


4U 1608-52 


619 


1.3 


0.35 


1 


14.5 


7.4 - 9.4 


4U 1608-52 


619 


1.4 


0.45 


1 


15.1 


6.9 - 8.9 


4U 1608-52 


619 


1.49 


0.4 


1 


16.5 


6.6 - 8.5 


slow rotator 




M 


{Rg/RY 


n, m 


R 


range of tqpo 


4U 1702-429 


330.55 


1.4 


0.4 


2 


28.0 


7.8 - 8.6 


4U 1702-429 


330.55 


1.0 


0.5 


2 


18.8 


10.0 - 11.0 


4U 1702-429 


330.55 


1.1 


0.4 


2 


23.0 


9.4 - 10.3 


4U 1702-429 


330.55 


0.83 


0.35 


2 


18.5 


11.6 - 12.7 


4U 1728-34 


364.23 


1.4 


0.4 


2 


22.5 


7.2 - 9.9 


4U 1728-34 


364.23 


1.2 


0.5 


2 


16.0 


8.1 - 11.1 


4U 1728-34 


364.23 


1.1 


0.4 


2 


16.5 


8.7 - 11.8 


4U 1728-34 


364.23 


1.1 


0.35 


2 


17 


8.7 - 11.8 


4U 1728-34 


364.23 


1.0 


0.4 


2 


17.5 


9.2 - 12.6 


estimated u s 


Sco X-1 


300.0 


1.4 


0.7 


2 


21.3 


7.4 - 8.8 


Sco X-1 


540.0 


1.4 


0.6 


1 


15.5 


7.0 - 8.0 


Sco X-1 


422.0 


0.9 


0.5 


1 


17.4 


9.9-11.5 


Sco X-1 


540.0 


1.2 


0.4 


1 


16.0 


7.9 - 9.1 


Sco X-1 


280.0 


0.8 


0.5 


2 


17.5 


11.4 - 13.3 


Sco X-1 


292.0 


0.81 


0.35 


2 


18.6 


11.3 - 13.2 



Note. — u s is given in unit of Hz, M in Mq, R in km, radial coordinate where QPO occurs, tqpo, in unit of GM/c 2 



Table 2 

Physical parameters of black holes 



black hole 


M 


estimated a 




vi 


tqpo 


Ar 


n, m 








theory / observation 


theory / observation 






GHO J 1655-40 


6-7 


0.737 - 0.778 


450/450 


300/300 


4.93 - 4.25 


1.71 - 1.23 


2 


GRO J 1655-40 


7.05 


0.95 


451.31/450 


299.04/300 


1.94 





1 


XTE J1550-564 


8-11 


0.682 - 0.768 


276/276 


184/184 


5.91 - 4.4 


2.44 - 1.33 


2 


XTE .11550-564 


11.5 


0.95 


276.67/276 


183.32/184 


1.94 





1 


H1743-322 


11.3 


0.74 


242.39/242 


165.89/166 


4.81 


1.6 


2 


H1743-322 


8-15 


0.644 - 0.82 


242/242 


166/166 


6.5 - 3.5 


2.86 - 0.7 


2 


H1743-322 


12.7 


0.95 


242.46/242 


166.05/166 


1.94 





1 


GRS 1915+105 


10 - 20 


0.606 - 0.797 


168/168 


113/113 


7.38 - 3.9 


3.58 - 0.98 


2 


GRS 1915+105 


18.4 


0.95 


167.35/168 


114.61/113 


1.94 





1 



Note. — are given in unit of Hz, M in Mq, tqpq and its distance from marginally stable orbit, Ar, are expressed in unit of GM/c 2 . 



O m - z 




Fig. 1. — Cartoon diagram describing coupling of various modes in an accretion disk and the corresponding nonlinear oscillator. The 
oscillators describing by the spring constant k\ and indicate respectively the coupling of spin frequency v* of the compact object with mass 
M to radial (x) epicyclic frequency and vertical (z) epicyclic frequency of a disk blob with mass m, where b\ and 62 represent corresponding 
damping factors respectively. 
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Fig. 2. — Variation of (a) radial (solid curve) and vertical (dotted curve) epicyclic frequencies in Hz as functions of radial coordinate in 
unit of GM/c 2 , (b) higher (solid curve), lower (dotted curve) QPO frequencies, their difference (dashed curve), u z — v r (long-dashed curve) 
in Hz as functions of radial coordinate in unit of GM/c 2 , for the parameter set given in Table 1 for KS 1731-260 when M = I.IMq and 
R = 16.4km. (c) Same as of (a), (d) same as of (b), but for 4U 1702-429 given in Table 1, when M = M and R = 18.8km. The vertical 
dot-dashed curve indicates the location where the resonance occurs to generate QPO. 
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Fig. 3. — Variation of (a) radial (solid curve) and vertical (dotted curve) epicyclic frequencies in Hz as functions of radial coordinate in unit 
of GM/c 2 , (b) higher (solid curve), lower (dotted curve) QPO frequencies, their difference (dashed curve), u z — v r (long-dashed curve) in Hz 
as functions of radial coordinate in unit of GM/c 2 , for the parameter set given in Tabic 2 for GRO J1655-40 when M = 6Mq, n,m = 2. (c) 
Same as of (a), (d) same as of (b), but for M = 7.05Mq, n,m = l. The vertical dot-dashed curve indicates the location where the resonance 
occurs to generate QPO. 
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Fig. 4. — Variation of the QPO frequency difference in a pair as a function of lower QPO frequency for (a) KS 1731-260, (b) 4U 1636-53, (c) 
4U 1608-52, (d) 4U 1702-429, (e) 4U 1728-34, (f) Sco X-l. Results for parameter sets given in Table 1 from top to bottom row for a particular 
neutron star correspond to the solid, dotted, dashed, long-dashed, dot-dashed (for 4U 1608-52, 4U 1728-34, Sco X-l only), dot-long-dashcd 
(for 4U 1608-52, Sco X-l only) lines. The dot-dashed and dot-long-dashcd lines in (c) and the dot-dashed line in (e) are to fit a part of the 
observed data points discarding the remaining ones. The triangles are observed data points along with error bars. 
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Fig. 5.— (a) Normalized Power-Spectrum for 4U 1636-53 showing 775Hz QPO, when M = 1.18M , R = 16.8km, (R G /R) 2 = 0.35, 
mechanical resistance r m = dv r /dr ~ 50scc — . (b) Variation of the normalized power in (a) as a function of radial coordinate in unit of 
GM/c 2 . (c) Variation of Q-factor as a function of mechanical resistance of the flow in unit of sec -1 . 



